Bulk Spin Polarization of Co( 1 _ x .)Fe x S2 
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We report on a new method to determine the degree of bulk spin polarization in single crystal 
Co( 1 _ ;E )Fe a ;S2 by modeling magnetic Compton scattering with ab initio calculations. Spin-dependent 
Compton profiles were measured for C0S2 and Coo.9Feo.1S2. The ab initio calculations were then 
refined by rigidly shifting the bands to provide the best fit between the calculated and experimental 
directional profiles for each sample. The bulk spin polarizations, P, corresponding to the spin- 
polarized density of states at the Fermi level, were then extracted from the refined calculations. 
The values were found to be P = —72 ± 6% and P = 18 ± 7% for C0S2 and Coo.9Feo.1S2 respec- 
tively. Furthermore, determinations of P weighted by the Fermi velocity (vf or v F ) were obtained, 
permitting a rigorous comparison with other experimental data and highlighting the experimental 
dependence of P on up. 

PACS numbers: 71.20.Be, 72.25.Ba, 75.50.Cc, 78.70.Ck 



Spintronic materials are the subject of considerable re- 
search and rapid technological development. In these 
materials and devices, the goal is to exploit the electron 
spin for applications such as data storage and read-heads. 
Some commercial devices exist, for example magnetic 
random access memory devices (MRAM) being based 
on the concept of spin-dependent tunneling magnetore- 
sistance [T]. Because of the potential for the develop- 
ment of novel devices using spin transport in semicon- 
ductors, so-called spin-injector materials (which would 
act as a source of highly spin-polarized electrons) cur- 
rently attract considerable interest. An obvious choice 
for a spin injector is a ferromagnetic system since it nat- 
urally has an imbalance of the electron spin population 
at the Fermi level and thus possesses a degree of spin 
polarization. Ideally, materials for applications would be 
fully spin-polarized 'half-metals' (HMs), where the den- 
sity of states (DOS) at the Fermi energy (Ep) is finite 
for one spin but zero for the other, such that carriers of 
only one spin exist at the Fermi level. However, ferro- 
magnets typically only possess partial spin polarization. 
Despite much theoretical and experimental effort, only a 
few promising candidates, such as Cr0 2 [2J[3], optimally 
doped Lao. 7 Sro.3Mn0 3 g] and Co 2 MnAl (1 _ :I . ) Sn a; [SJ 1] 
have been found to exhibit high polarizations. Recently, 
Co(i_ a; )Fe a ;S2 has garnered interest resulting from pre- 
dictions of half-metallicity and the potential to tune the 
polarization via Fe doping. The consensus is that doping 
simply adjusts the Fermi level p] [HI GO EH EE] , without al- 
tering the bandstructure. The possibility of altering the 
magnitude of the polarization makes this material ideal 



for fundamental research on spin polarization. 

The spin polarization can be defined simply in terms 
of the spin-dependent density of states (DOS), N^n, at 
the Fermi level. However, in order to facilitate compari- 
son with various experimental techniques, it is useful to 
expand this definition to weight iVj /j with the Fermi ve- 
locity, Vp^n. Hence, following Mazin [IS], P n can be 
defined such that, 

= N,v h - N ivh 

For n — 0, as quoted by most theoretical studies, Pq is 
solely defined by the DOS. Weighting with the Fermi ve- 
locity accounts for transport properties, with n — 1 in the 
ballistic or n — 2 in the diffusive regimes. Experimentally 
the direct measurement of the polarization has proven to 
be a challenge. In particular, except for the case of a HM 
where P n is 100%, the appropriate value of n is often am- 
biguous, and consequently comparison with theory is pos- 
sible, but notoriously difficult [H EH EH E3 ■ In point- 
contact Andreev reflection (PCAR) measurements, for 
example, the value of n is not the same for all materials 
|14j . Bulk techniques such as saturation magnetization 
or transport measurements can only give an indication of 
half-metallicity, and techniques which give access to the 
magnitude of P directly such as spin polarized photoe- 
mission or PCAR are heavily dependent on the surface 
quality of the sample. 

In this Letter, we demonstrate that a magnetic Comp- 
ton scattering study combined with ab initio electronic 
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structure calculations can be used to extract P n in the 
representative pyrite-type series Co(i„ K )Fe a; S2. Fine tun- 
ing the theory to experiment enables the spin polariza- 
tion to be determined for n — as well as the weighted 
values for n — 1, 2. 

Magnetic Compton scattering is a valuable tool for the 
investigation of spin-dependent effects as the incoming 
hard x-rays directly probe the twice-integrated momen- 
tum distribution of the unpaired electrons in the bulk 
|16j . A Compton profile (CP) represents the ID pro- 
jection of the electron momentum density. If the differ- 
ence between the CPs is measured for opposite applied 
magnetic field directions, the resulting difference profile 
contains only the spin-dependent part of the momentum 
distribution; this is referred to as the Magnetic Compton 
Profile (MCP). 

The directional MCPs can be directly compared to 
the spin-dependent momentum distributions which have 
been computed from the ab initio electronic structure 
calculation (from which the the polarization may also 
be calculated). It has been demonstrated previously 
[171 118j that the agreement with experiment may be re- 
fined by iterating a rigid shift of the bands (with respect 
to the Fermi level) contributing to the spin moment. In 
the case described here, this optimization is performed 
by minimizing the x 2 calculated simultaneously between 
the computed and measured MCPs along four crystal- 
lographic directions. The bulk polarization P, and the 
weighted polarization P n can then be calculated directly 
from the refined model bandstructure giving a direct ex- 
perimental insight into not only the polarization but also 
the effect of the Fermi velocity on P. The ability to ex- 
tract all possible permutations of the polarization allows 
simple comparison with results from other studies. 

In order to test our method, we calculated the spin 
polarization of Ni. For Ni, a transport measurement de- 
termined P2 to be 23 ± 3% experimentally [T5], whereas 
ab initio calculations predict P 2 = 0% [12] . For our re- 
finement, we used experimental MCPs for four crystallo- 
graphic directions from Dixon et al. |20j . The refinement 
increased the value of P2 from our ab initio value of 0% 
to 20 ± 2%, which agrees well with the previous experi- 
mental result. This simple test shows that it is a robust 
method for determining P n 

Bandstructure calculations predict C0S2 to be a highly 
spin-polarized ferromagnet, just missing half-metallic be- 
havior [21 . The Fermi level lies low in the conduction 
bands and doping with the isostructural semiconductor 
FeS2 results in a reduction in the occupation of one spin 
channel. It has been predicted that this would enable the 
polarization to be tuned to half metallicity [§] [TU] ITT] . 
Saturation magnetization and transport measurements 
indicate a high spin polarization hinting towards possi- 
ble half-metallicity at doping levels between 7 — 10%. 
The trend in the total DOS at Ep extracted from heat 
capacity and nuclear magnetic resonance follow this con- 




FIG. 1: (Color online) The DOS of C0S2, where the sul- 
phur p- and cobalt d-manifold are depicted for the majority 
and minority spin states. The inset shows the bandstructure 
of the e 9 -manifold around the Fermi level along a high sym- 
metry path in the cubic Brillouin Zone. Overall six bands 
composed of two in the majority (blue, solid lines) and four 
in the minority (red, dashed lines) channel cross the Fermi 
level. 



jecture [5] [22]. However, PCAR measurements of the 
absolute value of the magnetization reveal magnitudes of 
up to \P\ = 64% in the pure (single crystal, [23]) and 
\P\ = 85% for x = 0.15 (polycrystalline, p]) in the 
doped system. Complementing these results, anisotropic 
magnetoresistance measurements show a change in sign 
of the polarization on doping, as predicted by ab initio 
calculations [H] [24] . Although the qualitative behavior 
of both experiment and theory are similar, the absolute 
values of P disagree [25 . This may be attributed in part 
to the sensitivity of PCAR to (such that n > 0), as 
opposed to the theoretical values which refer to Pq. 

For the current study, measurements were performed 
on C0S2 and Coo.9Feo.1S2 single crystals to investigate 
polarization effects at low doping. High quality sin- 
gle crystals were prepared by the chemical vapor trans- 
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FIG. 2: (Color online) The DOS in the e 9 -manifold of (a) 
C0S2 and (b) Coo.9Feo.1S2 is composed of six bands around 
the Fermi level. As in Fig. [I] the DOS of the majority and 
minority bands are blue and red, respectively. The dark and 
light colors signify band 1 and band 2 (or in the minority 
channel, the combination of bands 2-4), respectively. 
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FIG. 3: (Color online) Experimental and LMTO MCPs for (a) C0S2 and (b) Coo.9Feo.1S2 resolved along different crystallo- 
graphic directions. For clarity the profiles are offset by 0.2 and statistical errors are smaller than the symbol size. The insets 
of both graphs show the directional differences for each composition between the [100] and [110]. 



port method [5] . Particular care was taken to ensure 
stoichiometry and to avoid sulphur deficiency, which is 
known to affect the measured polarization [25] . The 
magnetic Compton scattering experiments were carried 
out at the beamline (BL08W) of SPring-8, Japan. All ex- 
perimental data were collected at a temperature of 10K. 
The MCPs were subsequently corrected for detector effi- 
ciency, absorption, the relativistic cross-section and mul- 
tiple scattering effects. The value of the saturation mag- 
netization measured with a SQUID at 10K was deter- 
mined to be 0.83 hb/Co and 1.04 /z^/Co for the pure 
and x = 0.1 samples. 

The electronic stucture of C0S2 was computed using 
the linear muffin-tin orbital method (LMTO) within the 
atomic sphere approximation including combined correc- 
tion terms [27] . The exchange correlation part of the 
potential was described within the local spin density ap- 
proximation (LSDA) and the effect of doping was incor- 
porated via the virtual crystal approximation (VCA). As 
observed in previous calculations, the moment at the ex- 
perimental lattice constant of 5.5 A is substantially un- 
derestimated [551 HH] and so the calculations presented 
here have been performed at 5.71 A, producing moments 
which are 0.85 /is/Co and 0.92 /is/Co for the pure and 
x = 0.1 systems respectively, in good agreement with the 
experimental values. 

The results of the calculation for C0S2 are depicted 
in Fig. [T] The spin-resolved DOS show the sulphur 
2p-states below the cobalt 3c? manifolds. These states, 
which are hybridized with sulphur 2p-states, are split by 
the crystal field into t2 g - and e ff -manifolds. The Fermi 
level lies low in the e g -manifold where two majority bands 
and four minority bands cross Ep (see inset of Fig. [T]). 
The partial DOS contributed by each band are depicted 



in Fig. |2ja.) . In the undoped compound, in both the 
majority and minority spin channel, band 1 holds the 
largest DOS and is therefore expected to have the biggest 
impact on P n . The calculation predicts a polarization 
Po = —58%. The bandstructure and partial DOS for 
x = 0.1 are broadly similar to the pure case and lead 
to Pq = 40% (see Fig. [2]d). Our calculations predict 
half-metallicity to occur for x > 0.3. Overall these re- 
sults are consistent with other theoretical band predic- 
tions [241 Ell EH [30]. 

The experimental and calculated MCPs, resolved along 
different crystallographic directions, are shown in Fig. 
ga) and [gb) for CoS 2 and Coo. 9 Feo.iS 2 . The ab ini- 
tio theoretical results are in excellent agreement with the 
experimental MCPs. The characteristic structures in the 
directional differences of the magnetic Compton profiles 
serve as a rigorous check on the calculation, and these are 
shown for the pairs of directions [100] - [110] in the insets 
of Fig. [3] In order to optimize the agreement between 
theory and experiment the energy bands were refined, as 
described earlier, to provide the best fit to the experi- 
mental MCPs. Small energy shifts of the bands crossing 
Ep change the detailed shape of the resultant MCP, the 
band positions were refined to give the best fit between 
the theoretical and experimental MCPs, simultaneously 
for four crystallographic directions for C0S2 and three for 
Coo.9Feo.1S2. The spin moments calculated from these 
fits agree well with the experimental values, providing 
an independent consistency check. 

Despite small energy shifts of the bands (20 — 60 meV) , 
the fitting procedure sees the x 2 parameter reducing by 
30 — 50% from its starting value to the minimum and has 
an appreciable impact on the polarization. In the un- 
doped case, the polarization decreases from the original 



4 




O , □ ,-0- fitted P 012 
A- PCAR data from [9] 



0.00 0.05 0.10 0.15 0.20 0.25 0.30 
Fe doping x 

FIG. 4: (Color online) Experimental PCAR (filled triangles) 
polarization data compared to Po,i,2 from the LMTO calcu- 
lation (lines) as a function of composition in 5% steps. The 
tuned calculations are also shown, where the bare P is also 
weighted by the Fermi velocity (open symbols). Note that the 
sign change in the PCAR data has been inferred from AMR 
measurements. 



value of P = -58% to P a = -72(±6)%. Here the er- 
ror (which is generally larger for more dispersive bands) 
reflects the uncertainty of the fit. This decrease can be 
mainly attributed to majority band 2 being shifted above 
Fip and the DOS of band 1 being decreased. For the 
x = 0.1 compound, the polarization decreases from the 
original value of P = 40% to P = 18% (±7), which is 
primarily due to the lowering of majority band 2 and an 
increased DOS for the minority band. 

The fitted bulk polarization P n obtained for C0S2 and 
Coo.9Feo.1S2 are depicted in Fig. [4]along with the results 
from the unfitted calculation for different compositions. 
The results of the original calculation for Pq at different 
x are in good agreement with previous theoretical studies 
|24j . However, the weighting with the Fermi velocity has 
a significant effect on the nominal value of the polariza- 
tion, particularly in the pure system where a weighting 
with vp results in a positive polarization. The evolution 
of the fitted P n is qualitatively consistent with indirect 
measurements such as saturation magnetization, but im- 
ply that the maximum P is not obtained at x = 0.1. 

Our results demonstrate the tunability of the different 
polarization functions across the series of Co(i_ a; )Fe a ;S2. 
The data confirm that Pq is negative in the pure and pos- 
itive for the doped compound. It is instructive to com- 
pare our values for P n with the PCAR results. In the 
case of the pure compound the value of -PpcAR = |64%| 
is between Pq and P\. However, for the x = 0.1 com- 
position -PpcAR is m close proximity to P2 found in our 
study. This could indicate a doping dependent crossover 
from the ballistic to the diffusive regime in PCAR mea- 
surements, which is consistent with the increase in resis- 
tivity observed for higher Fe-doping 8 . Alternatively, 



these discrepancies could be due to surface effects play- 
ing a role in the PCAR measurements as pointed out by 
Leighton et al. [7]. Although PCAR has a penetration 
depth corresponding to the coherence length of the super- 
conductor, the electrons crossing the interfacial barrier 
experience the surface potential. Studies on the surface 
state of C0S2 show that sulphur degradation could influ- 
ence polarization measurements |31[I32], However, so far 
the surface stoichiometry of the alloyed composites has 
not been investigated. 

In summary, we present a new method to determine 
the bulk spin polarization. Although harder to put into 
practice, this novel technique calculates P n from a re- 
fined model bandstructure that is fitted to experimental 
MCPs. Moreover, it is free from surface effects. The case 
of Ni shows that this approach is a rigorous method of 
determining the bulk P n . More importantly, the results 
on C0S2 and Coo.9Feo.1S2 illustrate the applicability of 
the approach in a more complex system, and highlight 
the degree to which the polarization is dependent on vp- 
We speculate that there is a crossover from the ballistic 
to the diffusive transport regimes on doping with Fe on 
the basis of the comparison with PCAR measurements. 
The new method can be used to determine P n in such 
a way that its dependence on vp is unambiguous, and 
our results for Coi_ a: Fe x S2 demonstrate the importance 
of this. 
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